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ABSTRACT 



Aims. The previously unidentified very high-energy (VHE; E > 100 GeV) y-ray source HESS J1303-631, discovered in 2004, is re-examined 
including new data from the H.E.S.S. Cherenkov telescope array in order to identify this object. Archival data from the XMM-Newton X-ray 
" . . satellite and from the PMN radio survey are also examined. 

^ Methods. Detailed morphological and spectral studies of VHE y-ray emission as well as of the XMM-Newton X-ray data are performed. Radio 
• data from the PMN survey are used as well to construct a leptonic model of the source. The y-ray and X-ray spectra and radio upper limit are used 
to construct a one zone leptonic model of the spectral energy distribution (SED). 
?-H Results. Significant energy-dependent morphology of the y-ray source is detected with high-energy emission {E > 10 TeV) positionally coincident 
with the pulsar PSR J1301-6305 and lower energy emission (E < 2TeV) extending ~ 0.4° to the South-East of the pulsar. The spectrum of 
the VHE source can be described with a power-law with an exponential cut-off A^o = (5-6 ± 0.5) x 10"'^ TeV"' cm"^ s"', T = 1.5 ± 0.2) and 
Scut = (7.7 ± 2.2) TeV. The PWN is also detected in X-rays, extending ~ 2 - 3' from the pulsar position towards the center of the y-ray emission 
region. A potential radio counterpart from the PMN survey is also discussed, showing a hint for a counterpart at the edge of the X-ray PWN trail 
and is taken as an upper limit in the SED. The extended X-ray PWN has an unabsorbed flux of F2_iokev ~ l-6!^04 x 10"'^ ergcm"^ s"' and is 
detected at a significance of 6.5 cr. The SED is well described by a one zone leptonic scenario which, with its associated caveats, predicts a very 
low average magnetic field for this source. 

Conclusions. Significant energy-dependent morphology of this source, as well as the identification of an associated X-ray PWN from XMM- 
Newton observations enable identification of the VHE source as an evolved PWN associated to the pulsar PSR J1301-6305. This identification 
is supported by the one zone leptonic model, which suggests that the energetics of the y-ray and X-ray radiation are such that they may have a 
similar origin in the pulsar nebula. However, the large discrepancy in emission region sizes and the low level of synchrotron radiation suggest a 
multi-population leptonic nature. The low impUed magnetic field suggests that the PWN has undergone significant expansion. This would explain 
the low level of synchrotron radiation and the difficulty in detecting counterparts at lower energies, the reason this source was originally classified 
as a "dark" VHE y-ray source. 



Key words. Gamma-rays: observations - Pulsars: individual: PSR J1301-6305 - ISM: individual objects: HESS J1303-631 
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1. Introduction 

In recent years, nearly a hundred VHE -y-ray sources 
have been discovered by various experiments, including 
many different types of sources. Generally, sources from 
these different classes also exhibit radio and X-ray radi- 
ation, however, the discovery of TeV J2032-I-4130 by the 
HEGRA collaboration in 2002 ( |Aharonian et al.||2002l ) 
lead to a new class of extended Galactic VHE y-ray 
sources without obvious counterparts at other wavelengths. 
HESS J1303-631, serendipitously discovered by H.E.S.S. 
(High Energy Stereoscopic System) during an observation 



campaign of the pulsar PSR B1259-63 (Aharonian et al. 



2005a I, in observations taken between January and June 



of 2004 dAharonian et al.„2005b) . HESS J1303-631 was 
the first so-called "dark source" discovered by H.E.S.S. 
More of these sources were discovered by the H.E.S.S. col- 



laboration in the following years (Aharonian et al. 2008 



TiboUa et al. 2009) . Identifying and understanding this new 
class of sources has become an important task for modern 
y-ray astronomy. 

A growing number of extended VHE y-ray sources, 
without (or with significantly fainter or less extended) 
X-ray/radio counterparts, appear to be associated with 
energetic pulsars in the Galactic Plane. Some re- 
cent examples of this class of objects include HESS 
J1825-137 dAharonian etal.|2006aD and HESS J1356-645 
dAbramowski et al.l 120111). These associations are be- 



lieved to represent pulsar wind nebulae (PWNe), which are 
thought to be powered by a relativistic lepton-dominated 
particle outflow from a central pulsar. Many of these PWN 
associations form what are known as Offset PWNe where 
the pulsar is located at or near the edge of the y-ray and 
X-ray emission regions. These configurations may form in 
two ways. First, a high spatial velocity pulsar, possibly su- 
personic (in which case a bow shock nebula may form), 
leaves behind a "trail" of high-energy electrons in the am- 
bient medium. Alternatively, an offset PWN may form if 
its expansion is blocked on one side by the reverse shock 
of the supernova remnant (SNR) in which the pulsar was 
born. Due to inhomogeneous densities in the ISM, the ex- 
pansion of the supernova remnant may proceed asymmet- 
rically, or the motion of the pulsar may place it near the 
edge of the SNR and the expanding PWN may then be dis- 
rupted asymmetrically by the reverse shock of the SNR, a 
scenario known as a Crushed PWN ( |Blondin et al.|200T) ). 

At the time of discovery, HESS J1303-631 was found 
to have a large intrinsic Gaussian extent of ~ 0.16°, as- 
suming a 2-dimensional symmetric Gaussian distribution, 
and a flux of ~ 17% of the Crab flux above 380 GeV 
Originally, the source had no known extended counter- 
parts at other wavelengths and was, therefore, classified 
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as a dark source. As is the case with many such dark 
sources, HESS J1303-631 is found to have a pulsar lying 
near the edge of the emission region with a high enough 
spin-down luminosity to account for the y-ray emission. 
PSR J1301-6305 is located at the northwestern edge of 
the emission region of the H.E.S.S. source and, with a spin- 
down luminosity E = 1.7 x 10^^ erg s~^, is the most pow- 



erful pulsar within 6° of the H.E.S.S. source (Manchester] 
et al. ( 2005[ ), see Table[T]for a list of known pulsars within 



0.5° of HESS J1303-631). This pulsar is young, with a 
characteristic age of Tc = 1 1 kyr, and a rotation period of 
184 ms. 

Originally, the distance to PSR J1301-6305 was esti- 



mated to be 15.8 kpc, based on Taylor & Cordes (19931, 



which would have required a rather high y-ray conver- 
sion efficiency of 37% in the 0.3 to 10 TeV range. Using a 
newer model of the Galactic electron distribution, NE2001 



(Cordes & Lazio 2002 1, however, yields a much closer 



distance of 6.6 kpc. Using this updated distance estimate, 
the reported VHE y-ray spectrum at the time of discov- 
ery yields an integrated flux in the 1 to 30 TeV band of 
= l.2x 10~'^ ergcm"^ s~' or 3.7% of the current spin- 
down luminosity of this pulsar, (f 6.6 = E/47t{6.6 kpc)^ = 
3.26 X lO'^'^ergcm'^s"^), a y-ray conversion efficiency 
which is comparable to other VHE PWNe (typically 0-7%, 



see e.g. |Mattana et al.| ( |2009| )). 

A 5 ksec Chandra X-ray observation, partially cover- 



ing the VHE peak emission region (Mukherjee & Halpem 



2005 1, revealed several point sources within the field of 



view, but no extended emission corresponding to the y-ray 
emission region was found, and none of the radio pulsars in 
the field of view of the Chandra observation were detected. 
The possibility of an annihilating clump of dark matter as 
the origin of the y-ray signal was explored by [Ripken et al. 



(2008 1. Such a model could explain the lack of detection of 
lower energy counterparts. However, it was found that the 
spectrum obtained for this source would require an unrea- 
sonably high mass for the candidate dark matter particles 
(~ 40 TeV). Also, as mentioned in that study, the inferred 
lateral density distribution does not support a dark matter 
scenario. [Ripken et al. (2008 1, therefore, concluded this to 
be an unlikely candidate for the explanation of the VHE 
source. 

To build a complete picture of the y-ray emission pro- 
cess in this source, data from recent re-observations of 
HESS J1303-631 with the H.E.S.S. telescope array were 
analysed, enabling studies of energy-dependent morphol- 
ogy. Also, follow-up observations by the XMM-Newton X- 
ray satellite, performed in 2005, showing a detection of a 
compact source slightly offset from the pulsar position and 
a significantly extended PWN, are presented. 

In Section 2, the H.E.S.S. instrument, data and anal- 
ysis methods are discussed as well as the light curve. 
Section 3 describes the studies of energy-dependent mor- 
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phology of the y-ray source, followed by a discussion of 
the spectrum of the source in Section 4. Section 5 presents 
the results of the XMM-Newton X-ray follow-up obser- 
vations, showing an X-ray PWN associated with the pul- 
sar PSR J1301-6305. Finally, in Section 6, the implica- 
tions of the analysis are discussed and the case is made 
for an association of HESS J 1303 -631 with the pulsar 
PSRJ1301-6305. 



2. H.E.S.S. observations and analysis 

2.1. The H.E.S.S. instrument 

H.E.S.S. is an array of four imaging atmospheric 
Cherenkov telescopes located in the Khomas Highland 
of Namibia (23°16'18"S, 16°30'00"E) at an altitude 
of 1800 m above sea-level. The telescopes image the 
Cherenkov light emitted by charged particles in the exten- 
sive air shower created when a y-ray is absorbed in the at- 
mosphere. They are optimized for detection of VHE y-ray 
initiated showers in the energy range of hundreds of GeV 
to tens of TeV by Each telescope has a 107 m^ tessellated 
mirror surface and is equipped with a 960 photomultiplier 
tube camera with a field of view (FoV) diameter of ~ 5° 
( |Beml6hr et"ar][2003| [Comils et al."2003|). The telescopes 
are triggered in coincidence mode (Funket al.|2004,) assur- 
ing that an event is always recorded by at least two of the 
four telescopes allowing stereoscopic reconstruction of the 
showers. More information about H.E.S.S. can be found in 



Hinton (20041. 



2.2. Data and analysis tectiniques 

HESS J1303-631 was originally discovered during an ob- 
servation campaign for PSR B 1259-63. Follow-up obser- 
vations of the two sources between 2004 and 2008 led to a 
total dataset of 108.3 hours of live time, using only obser- 
vations which passed standard H.E.S.S. data quality selec- 
tion which rejects observations taken during periods of bad 
weather or with instrumental irregularities. The data were 
taken in wobble mode at an average zenith angle of 43.8°, 
with an average offset of 0.8° from the position reported in 
the discovery paper ( [Aharonian et al.|2005b 1. 

The data were analyzed using H.E.S.S. standard Hillas 



reconstruction (Aharonian et al. 2006b I. Cuts were ap- 
plied to the shower image parameters to minimize back- 
ground, primarily due to cosmic -ray protons. For spectrum 



extraction, standard cuts (also defined in Aharonian et al. 



(2006b)), were used together with the Reflected-Region 
Background method ( Berge et al.||2007 ) to subtract resid- 
ual cosmic ray background, which resulted in an average 
energy threshold of ~ 720 GeV. The resulting excess for 
this analysis was found to be 12085 photons for a detec- 
tion significance of 33 cr. Some of these early observations 
were made with telescope pointings coincident with the 
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Fig. 1. The HESS J1303-631 VHE y-ray excess map, produced 
using hard cuts and the ring background method, was smoothed 
with a Gaussian kernel with cr = 0.05°. Coordinates are J2000.0. 
The high spin-down power pulsar, PSR J 1301 -6305, is indicated 
with a green star to the North-West of HESS J 1303-631. The 
point source, associated to the pulsar PSR B 1259-63, is seen 
in the bottom of the FoV. The size of the H.E.S.S. PSF, also 
smoothed with a Gaussian kernel with cr - 0.05°, is shown in 
the white box to the lower left. The blue/red transition occurs at 
a detection significance of ~ 5cr. 



HESS J1303-631 emission region, rendering them unsuit- 
able for spectral analysis since placement of reflected re- 
gions for background estimation is not possible. For the 
morphology studies, hard cuts were applied to further re- 
duce background contamination and improve image recon- 
struction, and hence the point spread function (PSF) of the 
instrument, at the expense of a higher energy threshold, 
together with the Ring Background method, resulting in 
an average energy threshold of ~ 840 GeV. Cross-checks 



were performed using a multi-variate analysis ( Ohm et al. 



20091, with background suppression based on boosted- 
decision trees, leading to compatible results. 

2.3. VHE y-ray map and iiglit curve 

The VHE y-ray excess map (Fig. [TJ of the 
HESS J1303-631 FoV shows extended emission to 
the South-East of PSR J1301-6305. A fit of a two- 
dimensional asymmetric Gaussian function to the excess 
resulted in a best-fit position of a = 13'^02"^48^ ± 3^^^^, 
6 = -63°10'39" ± 24^;^t (J2000.0), with major/minor 
axis Gaussian widths of cr^ ^ 0.194° ± 0.008° and 
CTy = 0.145° ± 0.006°, with a position angle (counter 
clockwise from north) of ^ = 147° ± 6°. The x^lNDF 
of the fit was 390 / 391. The exposure gradient over 
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the source extension was found to be small and have a 
negligible affect on the resulting source position. The fitted 
position is consistent with the one quoted in the original 



discovery paper ( jAharonian et al.| |2005bj ), but slightly 
shifted towards the pulsar position due to the (compared 
to the discovery paper) higher energy threshold of the 
hard cuts used and the presence of energy-dependent 
morphology (see Section 3). 

The nightly flux was determined using a flux extrac- 
tion region of radius 0.6° to ensure full enclosure of the 
source, around the best fit position given above assum- 
ing a power-law spectrum with an index of 1.5. Studies 
were performed to account for influences from the nearby 
VHE source PSR B 1259-63. As expected for an extended 
source, with an estimated diameter of 40 pc at a distance of 
6.6 kpc, the nightly flux is consistent with constant emis- 
sion, with /NDF = 77/69, verifying the stability of the 
H.E.S.S. instrument over the period of data taking. 

3. Energy-dependent morphology 

To test for the presence of energy-dependent morphology 
in the VHE source, excess images were generated in the 
following energy bands: Ei = (0.84 - 2)TeV, E2 = (2 - 
10)TeV and E3 > lOTeV (Fig. |2} left, top to bottom). 
The radial acceptance of the FoV was determined from the 
data, thus naturally accounting for the energy dependence. 
Slices were made on the uncorrelated excess images hav- 
ing dimensions of 1 .0° x 0. 1 ° and centered at the best fit po- 
sition of the VHE excess. The orientation is chosen along 
the fitted position angle (see Sec. |2.2[). A Gaussian func- 



tion was then fit to each slice as shown in Fig. [2] (right). 
The intrinsic source width was obtained by fitting the con- 
volution of a Gaussian with the energy dependent H.E.S.S. 
PSE 

The resulting parameters of the PSF convolved 
Gaussian fits, mean c and intrinsic Gaussian width Wint, for 
each energy band (Table [2]l were then plotted as a function 
of energy (Fig. [3]). A fit of a constant to these parameters 
yielded very bad quality fits, which establishes the exis- 
tence of strong energy-dependent morphology. This mor- 
phology implies a spectral steepening in y-rays away from 
the pulsar, a physical property predicted to be present in 
evolved PWNe. Fitting a linear function yielded much bet- 
ter quality fits (Tab. |3]l and a model of the morphology pa- 
rameterized by a projected center of emission, c(E), cal- 
culated with respect to the pulsar position, and an intrin- 
sic source Gaussian width, Wint(£'), which is calculated by 
taking into account the (energy-dependent) finite angular 
resolution of the instrument: 



c = (0.157 ±0.012)° 
Wint = (0.215 ±0.012)° 



(0.006 ± 0.002)° X ^Tev 
(0.009 ± 0.002)° X ^Tev 



Tab le 1. All know n pulsars within 0.5° of HESS J1303-631 
( Manchester et al.l|2005^ . ^loiev is the distance from the given 
pulsar to the E > lOTeV peak position. 



Pulsar 

PSR J1301-6305 
PSR J1301-6310 
PSR J1305-6256 
PSR J 1302- 63 13 
PSR J1303-6305 



£/10^° erg/s 

1,700,000 

6,800 

760 

270 

7 



<5ioTcV [arc min] 
JTT 
6.6 
24.2 
7.2 
5.6 




1 Toen aw"! 3''02"i s'oo"' 

Right Ascension J2000.0 



/ 



1 Toen aw^ 3''02"i s^oo" 

Riglit Ascension J2000.0 




Distance trom puisar [°] 




Distance from puisar [°] 




1 3''06"1 3*^4"! 3''02"1 a'DO" 
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Distance from puisar [°] 



Fig. 2. Left: uncorrolated excess images of the HESS J 1303 -631 
region in the energy bands Ei = (0.84 - 2) TeV, E2 = (2 - 10) TeV 
and E3 > 10 TeV (from top to bottom). Coordinates are J2000.0. 
All images were smoothed with a Gaussian kernel of width 
0.05°. Slices are indicated by the rectangles, taken in the direc- 
tion of the semi-major axis of the fitted asymmetric Gaussian 
function. Right: the slices on the uncorrelated excess images 
are then fitted with a Gaussian function. The pulsar position is 
marked by a green star in the sky maps and a dashed line in the 
profiles. The dashed curves show the energy-dependent PSF of 
the H.E.S.S. instrument. 

Table 2. Results of the Gaussian fit to the slices on the excess 
images in the energy bands Ei - (0.84 - 2) TeV, E2 -(2 - 10) TeV 
and E3 > 10 TeV. c is the mean of the Gaussian, Wi,^„ is the 
Gaussian width and Wint is the intrinsic Gaussian width of the 
source after correcting for the PSF, P is the p-value of the x' fit- 
Band C W\„,„ Wint P 



El 
El 
Ei 



-0.14° 
-0.13° 
-0.05° 



:0.01° 
:0.01° 

: 0.02° 



0.21° 
0.17° 
0.08° 



:0.01° 
:0.01° 
: 0.02° 



0.20° 
0.16° 
0.07° 



:0.01° 
:0.01° 
: 0.02° 



0.42 
0.14 
0.97 
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Fig.3. Left: HESS J1303-631 fitted Gaussian mean, c(E), mea- 
sured from the pulsar position, as a function of energy. Right: 
the PSF corrected intrinsic Gaussian extension (wint(£), blue) is 
overlaid with the fitted uncorrected excess Gaussian extension 
(wimg(£'), black dashed) as a function of energy. The points are 
placed at the average energy of the photons falling in the corre- 
sponding energy bin (indicated by the horizontal error bars, not 
used in fit). 




E [TeV] 



Fig. 4. HESS J 1303 -631 spectrum derived using an integration 
region of radius 0.6°. The spectrum is well fit with a power-law 
function with spectral index F - 1.5 + 0.2 with a fitted cutoff 
energy of ^cut = (7.7 + 2.2) TeV. The fit resulted in a x^/NDF 
of (20 / 8) corresponding to 1% p- value. A fit of a power-law 
spectrum with no cutoff, shown by the dashed line, resulted in a 
p-value of 7 X 10"**. The last spectrum point deviates ~ 2cr from 
the fitted curve. It was removed from the residuals plot for better 
visibility. 



4. Energy spectrum 

The spectrum was derived using the Reflected-Region 
background method with an integration region of radius 
0.6°, roughly three times the intrinsic Gaussian extent at 
low energies to avoid effects of energy-dependent mor- 
phology, centered at the fitted source position. The de- 
rived spectrum for the entire dataset, excluding observa- 
tions where the offset of the pointing position to the center 



of the source is less than 0.6° (reducing the total live time to 
70.3 hours), is shown in Fig.|4] The spectrum was fit with a 
power-law function, dN/dE = No{E/l TeV)"'^, with a re- 
sulting photon index of F = 2.44±0.03stat and a normaliza- 
tion constant Nq = (5.9 ± 0.3stat) x 10"^^ TeV~^ cm'^ s"^. 
This normalization is larger than that found in the orig- 
inal discovery paper due to a larger integration region. 
However, with the inclusion of the additional data taken 
since the source discovery, the p-value of a chi-squared 
minimization is rather poor (7 x 10~^). A chi-square fit to 
a power-law spectrum with a cut-off at the energy Ecuu 



dN 
dE 



= NoE-^e'^'^"", 



yielded a better p-value of 1%, with fitted parameters Nq = 
(5.6 ± 0.5stat) X 10-12 TeV-icm-2s-\ F - 1.5 ± O.l^tat 
and ^cut = (7.7 ± 2.2stat)TeV. This spectrum yields an 
integrated flux in the 1 - 30 TeV band of (2.5 + 0.1) x 
lO"'' erg cm-2 s"' or 7.7% of Fg.e- Monte-Carlo studies 
were preformed to test for a possible contribution from the 
source PSR B 1259-63 (spill over events), due to the po- 
sition and size of the integration region and the exclusion 
region for PSR B 1259-63. Effects from this source are es- 
timated to be about 2% on the integrated flux, smaller than 
statistical and systematic errors. 

5. XMM-Newton X-ray observations 

In a search for counterparts of the VHE y-ray source 
in the keV energy band, two XMM-Newton observa- 
tions, each about 30 ksec, were carried out on July 
12* and 14"^, 2005, in satelUte revolution number 1024 
(ObsID 0303440101, "Observation 1") and revolution 
1025 (ObsID 0302340101, "Observation 2") respectively. 
All three X-ray imaging CCD cameras (EPIC MOSl, 
M0S2, and pn) were operated in full-frame mode, with 
a medium filter to screen out optical and UV light, with 
the exception of the pn camera during the first observation, 
where the Large Window mode with the Thinl filter was 
used. 

5.1. Data analysis 

For the data analysis of these observations, the XMM- 
Newton Science Analysis Software (SAS), version 9.0, was 
used (http://xmm.esac.esa.int/sas/). Cleaning the data and 
removing periods of high background due to soft proton 
flares resulted in a combined data set of about 52 ksec ex- 
posure. For this analysis, the energy band 2-8 keV was 
used to optimize the signal-to-noise ratio, since few events 
are expected at lower energies due to high absorption. The 
SAS task emosaicproc was used to combine the observa- 
tions and perform source detection, resulting in the de- 
tection of 73 point sources within the combined field of 
view above the maximum likelihood threshold of 10. The 
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Fig. 5. Top: The 2 - 8keV XMM-Newton X-ray count map in 
the region of the pulsar smoothed by a Gaussian kernel of width 
cr - 40". The horizontal axis is Right Ascension and vertical 
axis is Declination in J2000.0 coordinates. The count map is not 
exposure corrected, thus the apparent enhanced emission near 
the center is only an artifact. Overlaid is the projection annu- 
lus, shown in cyan, used to determine the direction of the X-ray 
extension, with inner radius of 48" and outer radius of 120", 
centered on the pulsar position. The 8, 14 and 20 cr TeV signifi- 
cance contours are shown in white. 

The direction of the extension is found to be within 1 cr from the 
direction of the star formation region, IRAS 13010-6254, lo- 
cated at a position angle of 106.3°, indicated by a magenta line 
in the sky map and a dashed line in the projection. The cyan 
lines in the sky map show the 1 cr errors in the fitted direction 
of the extension. Bottom: The X-ray azimuthal projection in po- 
sition angle from the pulsar location. The projected on-counts 
were fitted with the sum of a Gaussian and a flat background 
giving a position angle of 101.5° ± 5.3° and a Gaussian width of 
30° ± 7°. The point source located at ~ 210° is an unidentified 
X-ray source. 



X-ray PWN associated to PSR J 130 1-6305 was also de- 
tected in this way, corresponding to the 2XMMi catalog 
source 2XMM J130145.7-630536 dWatson et all 1200911, 



with a flux of F2-i2keV = (7.7 ± 1.0) X lO'^'^ergcm 
but peaked 15" ± 1.6" to the East of the pulsar with an ex- 
tension of 6" at a maximum likelihood of 7.7 (sources with 
likelihood < 8 may be spurious). No emission correspond- 
ing to the full extension of the H.E.S.S. source was found. 



However, a small apparently extended asymmetric X-ray 
source, directly adjacent to PSR J1301-6305, is seen ex- 
tending roughly 2' from the pulsar position towards the 
center of the VHE y-ray emission region (Fig. [5]). A de- 
tailed analysis of this feature is presented in the following 
section. 

While Observation 1 has the pulsar position closer to 
on-axis than Observation 2, it is unfortunately not suited 
for studying the extended X-ray source since the extended 
region found in Observation 2 lies directly on/between 
the edges of the CCD chips in all three detectors in 
Observation 1, thereby obscuring the view of this feature. 
Therefore, only Observation 2 was used for further analy- 
sis. 

5.2. Extended X-ray PWN 

To determine the direction of the X-ray feature, possibly 
associated to PSR J1301-6305, an annular projection was 
taken around the pulsar position with an inner radius of 
48" and an outer radius of 120" (Fig.|5]top). The projected 
counts were fitted with the sum of a Gaussian and a flat 
background giving a position angle of 101.5° ± 5.3° and a 
Gaussian width of 30° ± 7° (Fig. [5] bottom). The statistics 
were too low to warrant individual examination of the three 
cameras. The direction of the extension as determined here 
was used for the orientation of the slice on the count map, 
as presented below. 

The direction of the X-ray extension is consistent to 
within 1 cr with the direction of the star formation region 
IRAS 13010-6254, on the opposite side but within the 8 cr 
significance contour of the VHE source, as indicated in 
Fig. [5] top. This potential birthplace for the pulsar is con- 
sidered in more detail in Sec.|7] 

In order to determine the extension of the extended X- 
ray feature, a slice on the count map containing the pulsar 
was taken (Fig. |6j top) in the direction determined by the 
azimuthal projection, with a slice width of 88" and a length 
of 238" (on slice). A background slice of the same size and 
orientation was chosen in a source free region at roughly 
equal offset to the center of the FOV as the on slice to en- 
sure equal exposure. The slices are completely contained 
within single chips in the MOSl and M0S2 cameras and 
extend ~ 40" over the edges of neighboring chips in the pn 
camera. Profiles of the on slice and background slice are 
shown in Fig. [6] (middle and bottom). 

A point source located just West of the pulsar, 
2XMM J130141.3-630535, is presumably unrelated to the 
pulsar due to its soft nature (the hardness ratio R2, com- 
paring the 1 - 2keV to 0.5 - 1 keV bands, is 0.12 ± 0.08 
compared to 0.76 ± 0.12 for the source associated to the 
pulsar). 

The slice (Fig. |6j middle) does not exhibit enough 
statistics to precisely determine the morphology of the X- 
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Fig. 6. Smoothed, exposure corrected XMM-Newton X-ray flux 
map (cm^^ s ') in the 2-8 keV energy band (top) The horizon- 
tal axis is Right Ascension and the vertical axis is Declination 
in J2000.0 coordinates. An extended X-ray source appears from 
the position of the pulsar, PSR J 1301 -6305 (green dot), and ex- 
tends roughly in the direction of the center of the H.E.S.S. source 
(white contours show 8, 14 and 20 cr y-ray significance). Chip 
edges are shown in black for all three detectors and XMM cata- 
log sources are shown in magenta. The green box shows the slice 
used to create the profile (middle) and the red box shows the slice 
used for background estimation (bottom). Both slices were taken 
from the un-exposure corrected count map. A presumably unre- 
lated soft point source, 2XMM J130141. 3-630535, is shown in 
the map as a cyan dot. The excess was simultaneously fit with 
two Gaussians (the "compact region" near the pulsar position 
plus the "diffuse region" left of the pulsar) a King profile for 
the unrelated soft source to the right of the pulsar position and a 
constant background from the bottom slice. The red dotted bins 
in the projections lie directly on a chip edge in the pn camera 
in both slices and are excluded from the analysis. The dashed 
horizontal Une indicates the fitted background level. 



ray extension, but the extension appears to consist of a 
more compact region near the pulsar position, referred to in 
this section as the "compact region", and extending ~ 40" 
to the left, corresponding to the 2XMMi catalog source 
2XMM J130145.7-630536 as well as a feature extending 

from 40" to 150", referred to here as the "diffuse" 

emission region. 

A simultaneous fit of the slices was performed, con- 
sisting of a fit to the unrelated soft point source to the 
West of the pulsar, a Gaussian to the "compact" region 
near the pulsar position, a larger Gaussian to the "diffuse" 
region extending to the East and a constant to the counts 
in the background slice. The fit resulted in a diffuse emis- 
sion centered at -104" ± 18" with a Gaussian width of 
cr = 66" ± 19" while the compact region was found to 
be centered at -10" ± 4" with a width of 16" ± 4". The 
unrelated point source was fitted with a King profile 



c 



(1) 



with Rq - 4.3" and a - 1.5, corresponding to the PSF of 
the XMM-Newton pn camera at 1.5 keV and at ~ 10' offset 
from the center of the field of view. For the other cameras, 
the PSF is slightly narrower than this. 

The total X-ray extension is found to extend roughly 
170" (diffuse center -i-lcr width) beyond the pulsar posi- 
tion, however, the tail of the extension may be cut short 
by the edge of the pn chip. However, taking an integra- 
tion region from the edge of the pn chip, to avoid effects 
of changing sensitivity across chips, to the pulsar position 
(avoiding the soft point source to the west), for a total inte- 
gration length of 145" gives total on-counts of 950, and to- 
tal background counts 689 with the on/off area ratio a = I 
for an excess of 261 corresponding to a detection signifi- 
cance of 6.5 cr. 

5.3. X-ray spectrum 

For the spectral extraction, placement of a ring for back- 
ground determination was not possible due to multiple 
nearby sources. A small elliptical region covering the ex- 
tension region was taken and an identical ellipse was used 
for background extraction (Fig. |7]l. The extraction regions 
are well contained within single chips for the MOS and pn 
cameras. The spectra were obtained for the three cameras 
indepedently and then fit simultaneously. 

The obtained spectrum was fit using the spectral fit- 
ting software XSPEC with an absorbed power-law model 
which yielded a column density Nh = x 10^^ cm~^, 

a photon index T = 2.0^Qy, and a flux normalization at 
1 keV of 6.2+^° X 10"^ keV"^ cm'^ s'^ The integrated un- 
absorbed flux in the 2 - lOkeV energy band was found to 
be 



Fi-wkev = l.elli X 10*"ergcm-^s 
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Fig. 7. Zoom of Fig. [e] to the region around PSR J1301-6305. 
The region used for extraction of the X-ray spectrum is marked 
by the green ellipse. The red ellipse shows the region used for 
background determination. The VHE 14 and 20 cr significance 
contours are shown in white. 



6. Radio observations 

The region of HESS J1303-631 was covered by a sur- 
vey of the southern sky by the Parkes, MIT and NRAO 
(PMN) radio telescopes at 4.85 GHz ( [Condon et al.|1993| ). 
Calibrated maps were obtained from the NASA Sky View 
online tool, shown in Fig. [8] There is a radio feature just 
East of the X-ray nebula and near the peak of the VHE 
source, the apparent position of which may be shifted 
slightly to the North-East due to a strong gradient in the 
FOV from the strong unidentified radio sources to the 
North-East. The feature is found to have a peak flux of 
0.03 Jy/beam. The flux resolution (rms) of the PMN sur- 
vey is 0.01 Jy/beam so that the significance of this feature is 
only 3 cr and is at the detection limit of the survey (and thus 
not reported in the catalog). Therefore, the flux is taken as 
an upper limit. The feature is consistent with the size of the 
PSF of the survey (7' FWHM) in the North-East to South- 
west direction, but may be slightly elongated in the North- 
west to South-East direction, roughly parallel to the X-ray 
extension. Since the feature is not significant, no definitive 
conclusions about its morphology can be made. 

Although it is unclear whether this radio feature does 
indeed represent a counterpart of the y-ray and X-ray 
sources, since this lies in a rather complicated region of 
the radio sky, the location is promising due to its similari- 
ties with other known PWNe having a radio peak just be- 
yond the X-ray nebula, as in, for example, PSR B 1929-1- 10 
( [Becker et al.|2006 1 and the much smaller scale example of 
the "Mouse" PWN ( [Gaensler et al.|2004[ ). Further observa- 
tions in radio may be able to determine more precisely the 
morphology and polarization of the feature in order to es- 
tabUsh an association. 
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Fig. 8. 4.85 GHz radio image from the PMN Survey (beam size 
FWHM 7') in the HESS J 1303-631 region. The horizontal axis 
is Right Ascension and the vertical axis is Declination in J2000.0 
coordinates and scale is Jy/beam. H.E.S.S. contours are shown in 
green, XMM-Newton X-ray contours are shown in black and the 
radio contours are shown in white. A radio feature peaks about 
3' East of the pulsar position, just beyond the extended XMM- 
Newton X-ray source and near the center of the H.E.S.S. y-ray 
source at a peak value of 0.03 Jy/beam. The apparent position of 
the radio feature may be slightly shifted to the North East due to 
a strong gradient in background from neighbouring sources. 



7. Discussion 

7.1. Multi-wavelength situation 

Having analysed the morphology and spectra in VHE y- 
rays. X-rays and radio data available for the region, it is 
now possible to consider HESS J1303-631 in a full multi- 
wavelength scenario. First, an energy mosaic of the VHE 
emission was created using the three smoothed excess im- 
ages from Fig. |2] These images were overlaid, as shown in 
Fig.|9j along with the contours of the extended X-ray PWN. 
This energy mosaic is rather reminiscent of the known off- 
set PWN association HESS J 1825- 137 ([Aharonian et al. 



2006a I where the low-energy VHE y-ray emission region 
is quite extended with the pulsar laying towards the edge of 
emission and with the higher energy emission more com- 
pact and found centered closer to the pulsar. 

Taking the spectra and fluxes obtained in previous sec- 
tions, it is now possible to consider the SED of the source 
in a PWN scenario. Although a time-dependent model, in- 
cluding the evolution of the lepton populations over time, 
would be required to accurately describe the emission ob- 
served in the various wave bands, for simplicity, and due 
to the limited number of multi-wavelength data available, 



a simple stationary "one zone" leptonic model ( ) Aharonian 
& Atoyan 1999) was used to fit the VHE y-ray and X- 
ray spectra as well as the single PMN upper limit in radio 
(Fig. [T0| ). The leptonic model assumes that the same elec- 
tron population, with an energy distribution in the form of 
a single power-law with an exponential cut-off, creates ra- 
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dio and X-ray emission via synchrotron emission as well 
as VHE y-rays via inverse Compton (IC) scattering on 
Cosmic Microwave Background photons. Inclusion of IC 
scattering on infrared and optical target photons (as ob- 
tained from GALPROP (Moskalenko et al.|[2002 ]), assum- 
ing a pulsar distance of 6.6 kpc) had a negligible effect on 
the model parameters. 

The fit of the radio upper limit, and the X-ray and y-ray 
fluxes with this model yielded an electron spectral index of 



The VHE y-ray morphology presented here favors 
the association of HESS J1303-631 with the high spin- 
down power pulsar PSR J1301-6305, which, on energetic 
grounds, is the only known association which can explain 
the TeV emission. Additionally, the detection of an ex- 
tended asymmetric X-ray nebula in combination with a 
hint of a radio counterpart of the PWN seen in PMN ob- 
servations, strengthens this association further. 



a cut-off energy of Ecvn - 
ization of = 3.7+^-^ x 10^ cm'^ 



3 1!4 TeV, a normal- j^. Distance to the source 



and an average mag- 
netic field of 1 .4^Q^ juG, which is similar to the inferred 
mean line-of-sight magnetic field strength of ~ 2yuG, as 
determined from the pulsar's rotation measure ( [Crawford 



& Tiffany] 2007 1, but larger than the magnetic field of 
0.17 luG predicted by the y-ray to X-ray luminosity scal- 



ing law given in |Aharonian & Atoyan| ( |1999j ). The p- value 
of the fit was 0.02 and the model predicts a total energy 
in electrons of ~ 2.3 x 10^^ erg. It is worth noting that the 
resulting model spectrum in the radio band is steeper than 
typically observed in PWNe since the X-ray spectral index 
is not constrained and the fit of the single electron popu- 
lation is dominated by the narrow peak in TeV energies. 
Since the fluxes at the various energies described by this 
model are extracted from regions of differing size, the fit- 
ted magnetic field represents only an average and should 
be interpreted with caution. 

The differing sizes of the y-ray and X-ray emission re- 
gions imply the existence of differing electron populations 
so that the entire PWN cannot be accurately modeled by 
a single population. The simple approach presented here, 
therefore, suffers from the caveat that a model with two 
electron populations could reproduce the observed spectra 
with a significantly different magnetic field than obtained 
with a one zone model. Incorporating a strong cutoff in an 
older electron population at high energies would suppress 
the X-ray synchrotron emission, even in the face of a much 
higher magnetic field, and still reproduce the VHE peak. 
Indeed, the higher energy synchrotron emitting electrons 
may have been effectively extinguished precisely because 
of the high magnetic field. The morphology in VHE y-rays 
shows no evidence of a distinct break in the populations of 
electrons caused by passage of an SNR shock, but rather 
appears to show a more continuous transition from lower 
to higher energies in VHE y-rays and on up to the high- 
est energy synchrotron X-ray emitting electrons closer to 
the pulsar. This would imply a continuous transition from 
older to younger electrons which may require not a two 
zone electron model, but a continuously changing popu- 
lation making modeling quite difficult. Due to scant spec- 
tra available at lower energies, the precise details of the 
electron populations cannot be distinguished, and this first 
order approximation model serves as a starting point for 
future studies and searches. 



As stated before, the distance of 6.6 kpc to 
PSR J 130 1-6305 is based on the dispersion measure 
using the model of electron distribution in the Galaxy, a 
method which is often considered unreliable. For example, 
the "Mouse" pulsar has been argued to be at a distance 
roughly twice that determined by its dispersion measure 
based on the ratio of neutral hydrogen atoms to free 
electrons along the line of sight of A^///DM = 85 which 
is much higher than the values seen for all other X-ray 
detected pulsars, for which typically we observe A^///DM 
w 5 - 10 dGaensIer et aL]|^b4l ). For PSR J1301-6305, 
using the column density obtained here, we have A^/z/DM 
~ 23, one of the highest known A^/z/DM ratios among 
PWNe, which could imply that the distance obtained from 
DM is an underestimate for this source as well. 

On the other hand, the star formation region 
IRAS 13010-6254, if considered as a potential birthplace 
of the pulsar, provides an alternative estimation of the dis- 
tance. The direction of the X-ray extension is found to 
be within 1 cr from the direction of the star formation re- 
gion, IRAS 13010-6254, located at a position angle of 
106.3°, the only other identified object within the VHE 
emission region besides lower energy pulsars and stars. 
IRAS 13010-6254 has been identified as a point source 



in GeV y-rays in Fermi observations by Wu et al. (2011 



which could indicate the presence of evolved massive stars, 
which are the progenitors to pulsars, with colliding winds, 
or the presence of an SNR within the star forming region, 
or potentially a counterpart to a part of the VHE source 
given the bulge in the VHE significance contours at the 
position of the star forming region (see Fig.[5j top). As the 
only star formation region within more than a degree of the 
pulsar, and in the absence of another plausible SNR associ- 
ation, this provides a plausible candidate for the birthplace 
of the pulsar. This could be similar to the case of the X-ray 
feature G359.95-0.04 which was identified by [Wang et al. 
( 2006 1 as a PWN and found extending in the direction of 
the young stellar complex IRS 13, which was suggested as 
a possible birthplace for the yet undetected pulsar. 

IRAS 13010-6254 has a kinematic velocity of Vslr = 
(33.4 ± 3.2) km/s ( [Bronfman et al.|[T996l ), which corre- 
sponds to a distance of ~ 12.4 - 12.9 kpc using the circular 
Galactic rotation model of [Brand & Blitz| ( [1993| | updated 
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with the Galactic structure parameters of Reid et al. (2009 1. 
This is nearly double the distance of 6.6 kpc based on DM, 
placing the source close to the edge of the Galaxy. This 
kinematic distance is corroborated by the measure of the 
column density from X-rays, which is larger than the total 
integrated Galac tic HI column density in th at direction of 
1.9 X 10^2 cm-2 ( [Dickey & Lockmanfl996] ). 



If the pulsar was bom in IRAS 13010-6254 then it 
would have traveled 0.28° or ~ 62 pc, implying a very 
high transverse velocity of ~ 5,000 km/s if the charac- 
teristic age of 1 1 kyr is taken as the true age. This age es- 
timate is, however, often considered to be unreliable and 
the true ages may differ by a factor of 2-3. The true age 
of the pulsar, assuming constant braking index, is given as 
( [Manchester & Taylor|1977] ) 



T = 



{n - 1)P 



1 



P 



n-l 



(2) 



The characteristic age Tc is calculated assuming a braking 
index of « = 3 (i.e. pure electro-magnetic braking) and that 
the birth period, Pq is much less than the current period, P. 
The braking index has only been reliably measured for a 
handful of young pulsars ( Alpar & Baykal||2006 ) and was 
found to be less than 3 in every case, with the extreme case 
of with n = 1.4 ± 0.2 ( Lyne et aL]|1996[ ) implying an age 
5 times greater than predicted by Tc if the assumption of 
Po « P still holds, or PSR JI734-3333 with an index of 
n - 0.9 ± 0.2 (Espinoza et al. 201 1 1 implying an age which 



cannot be calculated with the above formula. Dupays et al. 



(2010 1 calculated the effects of quantum vacuum friction 
on the spindown of pulsars and found a braking index de- 
creasing as 1 - (1 - no)e~^' with A = P/P -i- P/P for pe- 
riod P and predict the braking index of the Crab pulsar 
at birth of just below 3 and that it will faU to ~ 2 in the 



next 2 kyr from its current value of 2.5. Gaensler & Frail 
( 2000| ) proposed that "characteristic ages greatly underesti- 
mate the true ages of pulsars" based on proper motion mea- 



surements of PSR B 1757-24. On the other hand, Migliazzo 



et al. ( 2002[ ) used proper motion measurements to show 
that the pulsar PSR B 195 1-1-32 is likely ~ 40% younger 
than its characteristic age, implying a non-negligible birth 
period. 

The very large "darkness ratio" of y-ray to X-ray lumi- 
nosity (in the 1-30 TeV and 2-10 keV bands respectively) 
for this source of 156, makes this the darkest identified 
PWN to date (the darker HESS J1702-420, darkness ratio 



1,500, is now beUeved to be an SNR, |Giacani et al.| ( |20TT) ). 
This could imply a relatively old age for PSR J1301-6305. 
Mattana et al.| ( |2009| ) considered the darkness ratios of 



PWNe and PWNe candidates detected by H.E.S.S., and 
found a logarithmic scaling with the characteristic age of 
the assumed associated pulsar. This scaling law predicts 
a darkness ratio of 6 for a pulsar with characteristic age 
1 1 kyr, whereas the darkness ratio measured here would 



predict a characteristic age of 48 kyr, more than 4 times 
older. 

A factor of 3 to 5 in the age of PSR J130I-6305 
would bring the pulsar velocity down to ~ 1 , 600 
to 1,000 km/s. This high velocity is not unreasonable 
given the two component pulsar velocity model by 
|Arzoumanian et aL (2002 1 which predicts ~ 15% of 
all pulsars to have a space velocity greater than 1,000 
km/s, but would place PSR J 130 1-6305 among the 
fastest known pulsars, including the Guitar Nebula pul- 
sar (PSR B2224-F65), PSR B1953-F50, PSR B1800-21, 
PSR B1757-24 and PSR B1610-50 all beUeved to have a 



velocity of > 1,600 km/s (see Cordes & Chemoff (1998 1 ; 



Caraveo (1993 1). In the case that the pulsar is much older 



than 11 kyr, IC cooling may play an important role for 
the oldest electrons, i.e. those created nearest the place of 
birth, leading to strong energy-dependent morphology as 
observed here. 

Adopting the kinematic distance of IRAS 13010-6254, 
the integrated y-ray luminosity would represent about 
28% of the current spin-down luminosity of the pulsar 
(£i2.6 = £/47r(12.6kpc)2 = 8.9 x 10"" ergcm'^ s"^). 
This y-ray conversion efficiency is higher than for typical 
PWNe (< 10% for PWNe with known y-ray and X-ray 
luminosities, [Mattana et al. (2009 1). However, high y-ray 
conversion efficiency may not be unreasonable consider- 
ing the very high "darkness ratio" of for this source of 
156, which implies low synchrotron losses. Thus, an as- 
sociation of PSR J 130 1-6305 to the star forming region 
IRAS 13010-6254 cannot be ruled out, while the associ- 
ation is supported by the direction of the X-ray trail, the 
high ratio of A^///DM, the large absolute value of Nh con- 
sistent with the entire integrated galactic column density in 
that direction, the possible bump seen in the TeV emission 
at the location of IRAS 13010-6254, the recent detection 
of IRAS 13010-6254 by Fermi as well as the high ratio 
of y-ray to X-ray luminosity. This larger distance may also 
help explain the difficulty of detecting counterparts at other 
wavelengths. 

8. Conclusions 

PWNe now appear to constitute the largest class of 
Galactic VHE y-ray emitters. The first dark source, and 
considered "prototypical" dark source, TeV J2032-I-4130 



discovered by the HEGRA collaboration ( Aharonian et al. 



2002 1, was recently found to be "not-so-dark", after deeper 
X-ray and radio observations have revealed weak but sig- 
nificant counterparts (Butt et al. [2008 1 and the identifica- 
tion of a pulsar first in y-ray s ( Abdo et al.|2009) and then 
in radio ( [Camilo et al.||2009l ). The work presented here 
has successfully identified energy-dependent morphology 
in VHE y-rays as well as an X-ray PWN counterpart of 
HESS J1303-631, which now appears also to belong to 
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Table 3. Quality of fit of a constant vs. a line to the source in- 
trinsic Gaussian extension and mean, measured from the pulsar 
position, as a function of energy. The much improved p-values 
of the linear fits as compared to the constant fits indicate the 
presence of significant energy-dependent morphology. 
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Fig.9. Energy mosaic of HESS J1303-631. The horizontal axis 
is Right Ascension and the vertical axis is Declination in J2000.0 
coordinates. Red: E, = (0.84 - 2)TeV, Green: Ej =(2 - 10)TeV 
and Blue: E3 > lOTeV. The highest energy photons origi- 
nate nearest the pulsar, PSR J 130 1-6305 (marked by the green 
dot). The visible red corresponds roughly to the lO-cr signifi- 
cance contour of the entire source. XMM-Newton X-ray con- 
tours are shown in black. A potential birthplace for the pul- 
sar, IRAS 13010-6254, as indicated by the X-ray extension, is 
shown by a blue circle. 
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Fig. 10. Spectral Energy Distribution of HESS J1303-631 fitted 
with a simple stationary leptonic model. The required magnetic 
field is ~ 1.4 yuG. 



the "not-so-dark", or "synchrotron under-luminous" class 
of VHE y-ray sources having peak synchrotron energy 
fluxes that are much lower than the peak fluxes in the VHE 
regime. The observations presented here support the inter- 
pretation of this source as a large cloud of electrons, ac- 
celerated by the pulsar, which emit y-ray radiation through 
the IC mechanism. These electrons can have an IC emis- 
sion lifetime of the order of the pulsar age, and can, there- 
fore, reflect the total energy output of the pulsar since birth, 
while the X-ray part of the PWN, generated by higher en- 
ergy synchrotron emitting electrons with a much shorter in- 
teraction time, decreases rapidly in time and reflects only 



the more recent spin-down power of the pulsar (de Jager 



|et al.||2009] ). While an association of the pulsar with the 
star formation region IRAS 13010-6254 is far from clear, 
it has been shown that it at least cannot be ruled out. The 
larger distance obtained from IRAS 13010-6254 may ex- 
plain the very high absorption column density obtained 
from X-rays. The pulsar's true age, distance and origin re- 
main open questions, as well as the details of the underly- 
ing electron populations that are responsible for the multi- 
wavelength emission. Current efforts to extend the radio 
and X-ray measurements of this source will be crucial for 
a deeper understanding of the processes at play. 

Many other extended Galactic y-ray sources which 
were previously unidentified are also finding associations 
with pulsars and PWNe as this class of sources continues to 
expand. The results obtained here also support the hypoth- 
esis that this "not-so-dark" source may be understood in 
the context of very low magnetic field, possibly in combi- 
nation with a large distance to the source, causing relative 
extinction of the X-ray counterpart. 
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